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ABSTRACT 

The Via Lactea simulation of the cold dark matter halo of the Milky Way predicts the existence 
of many thousands of bound subhalos with masses above a few xlO 6 M Q , distributed approximately 
with equal mass per decade of mass. Here we show that: a) a similar steeply rising subhalo mass 
function is also present at redshift z = 0.5 in an elliptical-sized halo simulated with comparable 
resolution in a different cosmology. Compared to Via Lactea, this run produces nearly a factor of two 
more subhalos with large circular velocities; b) the fraction of Via Lactea mass brought in by subhalos 
that have a surviving bound remnant today within r2oo and with present-day peak circular velocity 
V max > 2 kms -1 (> 10 kms -1 ) is 45% (30%). Most of the Via Lactea mass is acquired in resolved 
discrete clumps, with no evidence for a significant smooth infall; c) because of tidal mass loss, the 
number of subhalos surviving today that reached a peak circular velocity of > 10 kms -1 throughout 
their lifetime exceeds half a thousand, five times larger than their present-day abundance and more 
than twenty times larger than the number of known satellites of the Milky Way; e) unless the circular 
velocity profiles of Galactic satellites peak at values significantly higher that expected from the stellar 
line-of-sight velocity dispersion, only about one in five subhalos with V max > 20 kms -1 today must be 
housing a luminous dwarf. Any mechanism that suppresses star formation in small subhalos must start 
acting early, at redshift z > 3; f) nearly 600 halos with masses greater than 10 7 M are found today in 
the "field" between 7-200 and 1.5r2oo, i.e. small dark matter clumps appear to be relatively inefficient 
at forming stars even well beyond the virial radius; g) the observed Milky Way satellites appear to 
follow the overall dark matter distribution of Via Lactea, while the largest simulated subhalos today 
are found preferentially at larger radii; h) subhalos have central densities that increase with V max and 
reach pdm = 0.1 — 0.3 Mq pc -3 , comparable to the central densities inferred in dwarf spheroidals with 
core radii > 250 pc. 

Subject headings: cosmology: theory - dark matter - galaxies: dwarfs - galaxies: formation - galaxies: 
halos - methods: numerical 
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1. INTRODUCTION 

In the standard cosmological paradigm of structure for- 
mation (ACDM), the universe is dominated by cold, col- 
lisionlcss dark matter, and objects like the halo of our 
Milky Way are assembled via the hierarchical merging 
and accretion of smaller progenitors. Subunits collapse 
at high redshift, and when they merge into larger hosts 
their high densities allow them to resist the strong tidal 
forces that act to destroy them. Indeed, numerical simu- 
lations have shown that CDM halos are lumpy, teeming 
with self-bound subhalos on all resolved mass scales. It is 
now well established that the predicted subhalo counts 
vastly excee d the number of observed satellites of the 
Milky Way (jMoore et alJll999al iKlvpin et al1ll999f ). a 
"missing satellite" or "substructure" problem that has 
been the subject of many recent studies. While some 
models have attempted to solve the apparent small-scale 
difficulties of CDM at a more f undamental level, i.e. by 
reduc i ng small-sca l e pow er (e.g. iKami onkows ki fc Liddld 
l2000t IColm et "all |2000D . astrophysical solutions that 
quench gas accretion and star formation in small ha- 
los have also been proposed as a plausible way out. 

The latter category includes cosmic reionization (e.g. 

Via Lactea , the highest- resolution simulation of 
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Despite a wealth of ideas, however, a full characteri- 
zation of the substructure problem is limited by three 
main uncertainties: 1) the old tally of Milky Way lumi- 
nous satellit es is very incomplete, and is being rapidly 
revised (e.g. iBelokurov et al.ll2007D by the discovery of a 
new, large population of ultra-faint dwarfs in the Sloan 
Digital Sky Survey (SDSS); 2) in cosmological V-body 
simulations, the ability of subhalos to survive the hierar- 
chical clustering process as substructure within the host 
is particularly sensitive to resolution issues, as subha- 
los with numerically softened central densities are more 
easily disrupted by ti dal forces (e.g. iMoore et al.l 1 19961 : 
iDiemand et al.ll2007af ): and 3) a detailed comparison be- 
tween theory and observation in terms of, e.g., the peak 
circular velocity (V max ) of subhalos is hampered by the 
fact that Vniax for the observed Milky Way satellites is 
poorly cons trained by stellar velo city dispersion measure- 
ments (e.g. IStrigari et "aHl2007ah . 
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10 6 Mp, (like Leo IV, Co ma Berenices, Canes Venatici II, 
see ISimon fc Gehal l2007t ) : such small dwarf galaxies are 
typically resolved in Via Lactea with > 50 dark matter 
particle s, i.e. with a particle mass that is two dex smaller 
than in iMoore et al.l (fl999bl ) . three dex smaller than in 
iKlvpin et all (119991). a nd more than one dex smaller than 
in the lDiemand et all (|2004f) and lGao et all (|2004f ) simu- 
lations of galaxy halos. Via Lactea follows the formation 
of a Milky Way-sized halo in a WMAP 3-year cosmology: 
within its r2oo (the radius enclosing an average density 
200 times the mean matter value), more than 2000 sub- 
halos can be identified at z = with M sub > 4 x 10 6 M , 
distributed approximately with equal mass per decade of 
mass. Of these, 40 are found within 50 kpc, a region that 
appeared practically smooth in previous simulations. In- 
terestingly, the number of dwarf satellites in the inner 
halo of the Milky Way has recently tripled. Six dwarfs 
are currently known within 50 kpc of the Galactic Center: 
the LMC, Sagittarius, Urs a Majo r II, Coma Berenices, 
Willman 1, and Segue 1 (iMateol 119981: iWillman et all 
120061 IZucker et all 120061 : iBelokurov et al.ll2007f) . for an 
estimated total number (after correcting for the SDSS 
sky coverage) of order 25. Note, however, that Willman 
1 and Segue 1 lie at the uncertain boundary between dark 
matter-dominated dwarfs and globular clusters. 

The apparent conflict between the Galaxy's relatively 
smooth stellar halo and the extremely clumpy dark mat- 
ter distribution of Via Lactea raises a number of impor- 
tant questions. First and foremost, is the amount of sub- 
structure found in Via Lactea typical of galaxy halos? If 
so, how efficiently does gas accretion, cooling, and star 
formation need to be suppressed as a function of subhalo 
mass to explain the observed abundance of dwarfs? Be- 
cause of tidal mass loss, many of the subhalos that have 
small masses and circular velocities at the present epoch 
were considerably more massive and could have formed 
stars in the past. Does a single gas removal/star forma- 
tion quenching mechanism determine the nature/fate of 
the subhalo population? What subset of subhalos has 
properties similar to the known stellar subsystems of the 
Milky Way? How does a dwarf's luminosity relate to 
the mass evolutionary history of its host subhalo? Is the 
star formation suppression mechanism a strong function 
of the environment, or are small dark matter clumps in- 
efficient at forming stars even beyond the virial radius of 
the Milky Way? We address some of these issues here. 
The paper is organized as follows. In § 2 we compare the 
amount and distribution of Via Lactea substructure with 
that measured in an elliptical-sized halo simulated with 
comparable resolution in a different cosmology. In § 3 we 
show that, because of late tidal mass loss, the number of 
subhalos that reached dwarf-galaxy sizes before accretion 
onto Via Lactea is much larger than their abundance to- 
day and the abundance of known Milky Way satellites. 
If substructure mass regulates star formation, then for 
a given mass threshold many more subhalos should have 
been able to build a sizeable stellar mass at some point in 
the past than indicated by their present-day abundance. 
Sections 4 and 5 discuss the main constraints on mecha- 
nisms for the suppression of star formation in small sub- 
halos. The subhalo radial distribution and central dark 
matter densities are discussed and compared to the ob- 
servations in §6. Finally, we present a brief summary 



and our conclusions in § 7. 

2. CDM SUBSTRUCTURE IN GALAXY HALOS 

The Via Lactea sim ulation was performed with the 
PKDGRAV tree-code (|Stadelll200l1 ) and employed mul- 
tiple mass particle grid in itial conditions gen erated with 
the GRAFICS2 package (|Bertschingerll2001f) . The high 
resolution region was sampled with 234 million particles 
of mass m p = 2.1 x 10 4 M© and evolved with a force res- 
olution of 90 pc starting from redshift 50. We adopted 
the best-fit cosmolo gical parameters fr om the WMAP 
3-year data release (jSpergel et al.ll2007l ): fl M = 0.238, 
Oa = 0.762, H = 73 kms^Mpc" 1 , n = 0.951, and 
o"8 = 0.74. The z = host halo mass within r2oo = 389 
kpc is Mhost = 1-77 x 10 12 M Q . Its circular velocity pro- 
file reaches a peak of Vhost = 181 kms^ 1 at n/host = 62 
kpc. Note that the size of the Galaxy halo is still quite 
uncertain: since the effects of disk formation on the total 
mass distribution are still poorly constrained, the possi- 
ble range of dark matter halo sizes that could form and 
host a 220 kms -1 di sk is quite large, w ith Vhost between 
150 and 220 kms -1 (jKlypin et al.ll2002f ). When compar- 
ing Via Lactea to the Milky Way, this uncertainty allows 
one to rescale velocities and distances by ±25%, and sub- 
halo abundances above a given V max by about a factor 
of 1.25 3 = 2. A comparable uncertainty i n the abun- 
danc es comes from halo-to-halo scatter (e.g. iReed et "all 
I2005D . Throughout this paper we will for simplicity as- 
sume that Via Lactea matches the properties of the Milky 
Way halo closely, but the substantial uncertainties men- 
tioned above should always be kept in mind. 

For comparison purposes, we have run another simu- 
lation, termed "le8EU" , that follows the formation of a 
more massive host halo with Mh os t = 1-34 x 10 13 M Q , 
f200 = 512 kpc, Vhost = 390 kms -1 and r V host = 116 
kpc. The high resolution region of this second run 
was sampled with 136,432,000 particles of mass m p = 
1.8 x 10 5 M© and evolved from redshift 59.5 down to 
z = 0.47 with a force resolution of 162 pc. For lc8Ell 
the best-fit cosmolo gical parameters fro m the WMAP 
1-year data release (jSpergel et al.l I2003T ) were adopted: 
Q M = 0.268, fl A = 0.732, H = 71 kms" 1 Mpc" 1 , 
n = 1, and erg = 0.9. 

In both runs subhalos were identified using the phase- 
space friends-q f -friend s algorithm 6DFOF described in 
iDiemand et al.l (|2006f ) . Around the centers of all groups 
with at least 16 members, spherical density profiles were 
constructed, and the resultin g circular velocity profiles 
were fi tted with the sum of a iNavarro. Frenk. fc White] 
(|1997f ) profile and a constant density background. Sub- 
halo tidal mas ses were assigned accor ding to the defi- 
nition given in IDiemand etaD (|2007al ). The wealth of 
substructure resolved in our elliptical-sized halo is com- 
parable to Via Lactea and is clearly seen in Figure [TJ 
The halo finder algorithm identifies comparable num- 
bers of dark matter clumps within the virial radius of 
Via Lactea at z = and that of le8Ell at z = 0.47. 
Figure [2] shows the differential and cumulative subhalo 
mass functions within r2oo for both runs. Substructure 
extends down to th e smallest mass scales reliably re- 
solved: as shown bv IDiemand et all (|2007af ). numerical 
resolution effects start to flatten the mass distribution of 
subhalos below about 200 particles, and the maximum 
circular velocity distribution below V max / Vhost ~ 0.02. 




Fig. 1. — Projected dark matter density-square map of our simulated elliptical-sized halo ("le8EU") at z = 0.47. The image covers an 
area of 980 X 980 physical kpc, and the projection goes through a 980 kpc-deep cuboid containing a total of 120 million particles and 
25,000 identified subhalos. The logarithmic color scale covers 24 decades in density-square. 



In the range 200m p < M su b < 0.01Mi los t, the best-fit 
slope of the differential distribution, dN/dM su b oc M^ b , 
is a = -1.86 ± 0.02 for le8Ell and a = -1.90 ± 0.02 
for Via Lactea. In the same mass range the cumulative 
mass function has slope —0.92 for lc8Ell and —0.97 for 
Via Lactea. Both simulations are therefore characterized 
by steeply rising subhalo counts that, in the case of Via 
Lactea, correspond to approximately equal mass per sub- 
structure mass decade: most subhalos are of low mass. 
Figure [3] (left panel) depicts the fraction of the host halo 
mass within a sphere of radius r2oo that is bound up 
in substructure less massive than M sub , / su b(< M sub ) 
as a function of M sub /A/ host . We measure a total mass 
fraction in substructure that is about 5% in Via Lactea 
and exceeds 9% in le8Ell: its radial distribution can be 
approximated as / su b(< r) oc r for 0.1 < ?*/r2oo < 1. 
Because of the steepness of the cumulative mass func- 
tion, these fractions appear to be converging rather 
slowly at the small- mass end: more than 1% of the 



host mass is found in clumps with Ms U b/Mi lost < 10~ 5 . 
The amount of massive substructure is expected to in- 
crease with host halo mass since more massive hosts form 
later and dynamic al friction and tidal-stripping have less 
time to operate dGao et al.l l2004t iZentner et all 120051 : 
van den Bosch et al.ll2005fl . For the same reason, par- 
ent halos of a given mass will have a larger abundance 
of subhalos at higher rcdshifts than their present-day 
counterparts. However, when comparing z = clus- 
ters with galaxy halos these effects arc smaller than 
the observed halo-to-halo scatter in sub halo abundance 
(jDiemand et al.l 120041 : iReed et al.l l2005h . Variation in 
the slope and normalization of the power spectrum can 
also affect the amount of subhalos significantly, espe- 
cially their circular velocity function (jZentner fc Bullockl 
120031) . 

It is conventional to discuss the substructure popula- 
tion of galaxy halos in terms of the circular velocity func- 
tion. Figure [3] (right panel) shows the cumulative max- 
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Fig. 3. — Left panel: Fraction of the host halo mass that is bound up in substructure less massive than M au \, as a function of -Msub/A^host- 
Lower solid curve: Via Lactea. Upper dashed curve: le8Ell. Right panel: Cumulative maximum circular velocity function, N(> V ma x), 
for all subhalos within r200i as a function of V max / Vhost- There are 7932 subhalos with V max /Vh os t > 0-01 m Via Lactea, versus 8005 in 
lc8Ell. Also plotted is the radial distribution, 7V(< r), of all mass-selected subhalos with M su vj/-Mh st > 10~ 6 versus r/r200- Solid curves: 
Via Lactea, 7V(< r20o) = 4021. Dashed curves: le8Ell, N(< r' 2 oo) = 4447. 



imum circular velocity (Vmax = max{y / G?M su b(< ?')/''}) 
functions for the entire subhalo population within the 
two hosts. These are well-fit by a power-law of the form 

N{> K^x) = 100y- a 3 Xjl0 (0.4 < V m ax,10 < 3.5) (1) 

for Via Lactea, and 

A^(> Knax) = 1500 V ma 2 x 8 10 (0.8 < Vnax.lO < 5.5) (2) 

for le8Ell, with fractional errors not exceeding 20% in 
both cases. (Here Vmax.io = V ma x/10 kms -1 .) When 
Knax is normalized to the maximum circular velocity 
of the host, Vhost, the two substructure velocity func- 
tions are very similar at the low end, and become more 
discrepant at larger Vm ax values. There are 110 sub- 
halos above V ma x = 10 kms -1 = 0.056 Vhost hi Via 
Lactea, compared to 170 above Vmax = 22.5 kms -1 = 
0.056 Vhost hi lc8Ell. As already discussed above, it is 



well known that large subhalos (roughly within a fac- 
tor of 30 of the maximum circular velocity of the host) 
constitute a transient population that continuously de- 
clines from mass-loss effects. The Via Lactea run has 
now shown that small-scale substructure is instead more 
persistent and indepen dent of the age and size of the host 
l|Diemand et al . 2007b), consistent with pr edictions from 
semi-analytic models (jZentner et al.l I2005T ) . This is also 
supported by the left panel of Figure [3] within the range 
10 -6 < Af su b/Mhost 10~ 2 , the measured substructure 
mass fractions / S ub(< Msub) in the two simulations differ 
by less than ~ 40%. While the similarity between the two 
objects is noteworthy, we caution that such close agree- 
ment could be just a coincidence, and substantial amount 
of halo-to-halo scatter may be present in the abundance 
of small-scale substructure. 4 The larger dumpiness of 

4 Note that the measured abundance of small-scale substructure 
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lc8Ell is dominated by its three surviving subhalos with 
M su b/-Mhost > 0.006, versus none this massive identified 
in Via Lactea. 

The right panel of Figure [3] also shows the radial distri- 
bution of subhalos as a function of r/r 2 oo- In the range 
0.1 < x = r/r2oo < 1 ; the abundance profiles of mass- 
selected subhalos are very similar and well fit by 



N(< x) = N{< r 200 ) 



12x J 



1 + llx 2 



(3) 



This distributi o n agr ees well with those p resented by 
iDiemand et~aT1 (|2004l ) and iGao et all ()2004l ). extending 
the validity of these previous results down to M su b > 
10 _6 Mi lost . Subhalos do not trace the matter distri- 
bution of the host: tidal disruption is most effective in 
the inner halo, which leads to an antibias in the abun- 
dance profile of substructure relative to the smooth back- 
ground, i.e. the radial distribution of subhalos above a 
given mass is much flatter than that of the dark mat- 
ter in the host halo. This antibias is less pronounced if 
subhalos are selected by today's peak circular velocity 
instead. The threshold V max /Vhost > 0.02 yields a sam- 
ple of 2016 and 2334 subhalos in Via Lactea and le8EU, 
respectively. In the range 0.01 < x < 1, the abundance 
profiles of these V max - selected subhalos is well fit by 



N(< x) = N(< r 200 ) 



20a; 2 



1 + 19a; 1 



(4) 



3. FOSSIL RECORDS AND MISSING SATELLITES 



The Via Lactea and le8Ell runs discussed above show 
that, in the standard CDM paradigm, galaxy halos are 
filled with tens of thousands low-mass subhalos, about an 
order of magnitude more than found in previous simula- 
tions. Subhalos lose mass from the outside in through 
tidal stripping, the largest average ma ss fraction be- 
ing los t at the first pericenter passage (|Diemand et al.l 
l2007b[ K Larger clumps retain less of their mass as they 
sink quickly towards the host center due to dynamical 
friction, while the retained mass fraction is larger for 
initially lighter subhalos that are less affected by dynam- 
ical friction and stay in place. As many as 97% of all 
well-resolved [> 5,000 particles, V ma x(z) > 10 kms" 1 ] 
subhalos identified in Via Lactea at z = 1 do not lose 
their identity and have a surviving z = remnant. This 
can be compared with the 47% that survive from z = 4 
and the 20% from z = 10. Note that not all the subunits 
without a remnant today are totally disrupted to become 
part of the "smooth" host halo, some simply merge into 
larger substructure that eventually survives down to the 
present epoch. 

It is interesting to know, at the high-resolution limit of 
our simulations, what fraction of the Via Lactea's host 
halo mass was accreted in identifiable clumps, and com- 
pare it with the mass brought in by subhalos that have 
a surviving bound remnant within r2oo at the present 
epoch. The latter is expected to be significantly smaller 
than the former, since many major mergers at early times 
leave no surviving z = subhalos. We have there- 
fore traced backward in time the progenitors of present- 
day remnants in each simulation snapshot, measured the 

may also be sensitive to our change from a WMAP year- 1 (le8Ell) 
to year-3 (Via Lactea) cosmology, especially to the lower normal- 
ization of the power-spectrum as ■ 



maximum mass they had at any point during their his- 
tory i.e. before they suffered tidal mass losses, summed 
up this mass over all subhalos within r2oo having peak 
circular velocity > Vmax, and divided the result by the 
mass of Via Lactea today. The total mass fraction, 
/acc(> V max ), brought in by surviving substructure is 
plotted in Figure 2] as a function of Vmax, and is found to 
exceed 45% for V max > 2 km s" 1 . There are (14, 109, 873) 
subhalos today with V max > (20, 10, 5) kms -1 , which 
contributed (19, 30, 39)% of the mass of Via Lactea. Only 
9-10% of the total material they brought in remains in 
self-bound identifiable clumps at the present-epoch, the 
rest having been removed by tidal effects to become part 
of the "smooth" halo component. Note that numerical 
effects lead to less cuspy subhalos that suffer from arti- 
ficially high-mass loss, especially in the inner ha lo where 
tidal forces are strong ijKazantzidis et al1l2004f ). Figure 
[4] also shows the total mass accreted in clumps, f acc (> 
V maXiP ), independently of whether these survive or not to 
the present day. This is plotted against the highest Vm ax 
reached by subunits throughout th eir lifetime, a quan - 
tity we refer to as 'Vnax.p' (see also lKuhlen et al.ll2007l ). 
Interestingly, more than 97% of the Via Lactea mass is 
acquired in resolved discrete clumps with no evidence for 
significant smooth accretion. 

The large tidal mass losses suffered by massive subha- 
los have obvious implications for astrophysical solutions 
to the satellite problem. The leading scenario for rec- 
onciling the mismatch between luminous and dark sub- 
structure in the Milky Way posits that stars are able to 
form continuously and efficiently only in halos above a 
certain mass threshold, a threshold determined by the 
inability of shallow potential wells to accrete intergalac- 
tic gas photoionized by the ultraviolet (UV) background 
and/or to retain gas blown out by supernova explosions. 
There are slightly in excess of 20 dwarf galaxies known 
to orbit the Milky Way. After correcting for the sky 
coverage of the fifth data release of the SDSS, i.e. af- 
ter weighing each of the new ultra-faint dwarfs by a 
factor ~ 6, t he total number of satellites is probably 
around 60-70 (|Simon fc Gehall2007t iKoposov et al.ll2007h 
down to a surface brightness limit of fxy = 30 mag 
arcsec -2 . A solution to the subst ructure problem sue 



gested bv lStoehr etafl (|2002f ) and lHavashi et all (|2003l ) 
would be to place the luminous dwarfs in the most mas- 
sive subhalos at the present epoch, suitably adjusting 
the mass (or peak circular velocity) threshold for effi- 
cient galaxy formation. There are 65 subhalos today 
in Via Lactea (within r2oo) with M su b > 1.4 x 10 s M 
or Vmax > 12 kms" 1 : the paucity of Milky Way satel- 
lites would then be simply explained by the inhibition 
of star formation below the above values of M su t> and 
Vmax- Baryonic material falling in the dark matter po- 
tential well will be shock heated to the effective virial 
temperature of the host, compressed, and eventually cool 
and fragment into stars. A circular velocity threshold of 
12 kms" 1 corresponds to a virial temperature 



2k B 



< 8700 fxK, 



(5) 



where \i is the mean molecular weight. Figure [5] shows 
the cumulative number of Via Lactea subhalos as well 
as all Milky Way satellite galaxies within 420 kpc (the 
distance of Leo T, llrwin "eTaDl200l as a function of cir- 
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Fig. 4. — Left panel: Via Lactea accreted mass fraction / acc (> Vm ax ) (total mass brought in by all substructure surviving today within 
f*200 above V max , normalized to the mass of the host), as a function of today's subhalo circular velocity (upper solid curve). For comparison, 
we also plot the mass / aC c(> V ma x,p) brought in by all resolved subunits as a function of the highest maximum circular velocity reached 
throughout their lifetime, V max ,p (upper short-dashed curve). The long-dashed lower curve shows the mass fraction / su b(> V ma x) that 
remains in self-bound identifiable clumps at the present-epoch. Right panel: Cumulative circular maximum velocity function of all subhalos 
today within 50, 100, and 389 kpc (solid lines, from bottom to top). Dashed lines: Same for a sample of subhalos having a surviving bound 
remnant today and selected by their V ma x,p- 



cular velocity. The current available data are summa- 
rized in Table 1. The data point s in the figur e include 
all the previously known dwarfs (|Mateol[i998l ) plus the 
new circular velocit y estimates of the ult ra-faint Milky 
Way satellites from ISimon fc Gehal d2007h. plus Boote s 
(|Munoz et al.ll2006t ) and Willman 1 ([Martin et al.ll2007f ). 
They have been plotted assuming a maximum circular 
velocity of V max = \/3a ([Klypin et al.lll999l ). where a is 
the measured stellar linc-of-sight velocity dispersion, i.e. 
assuming a stellar spherical density profile oc r~ 3 in a 
singular isothermal potential. Note that the assumption 
of a constant multiplicative factor between V max and a is 
merely the simplest thing to do, and is not likely to hold 
on a case-by-case basis. Detailed modeling of the radial 
velocity dispersion profile, allowing for variations in the 
DM mass distribution and the stellar velocity anisotropy 
([Strigari et al.ll2"o"07af ). would be preferrable, but is cur- 
rently only available for a subset of all known dwarfs. 

If the stellar systems deeply embedded in dwarf 
sphcroidals remain largely unaffected by tidal stripping 
(this is clearly not the case, e.g., for Ursa Major II and 
Sagittarius), then the mass removal of large fractions of 
their original halo mass by tidal effects may make so- 
lutions in which luminosity tracks current subhalo mass 
somewhat misleading. Our simulations show and quan- 
tify better than before that many of the dark matter 
clumps that have small masses and circular velocities at 
the present epoch were considerably mo re massive and 
shoul d have formed stars in the past (e.g. iKravtsov et al.l 
12004 ) . We illustrate this point in Figure [4] (right panel) , 
which shows the cumulative circular maximum velocity 
function of substructure within 50, 100, and 389 kpc. 
Also plotted, for comparison, is the abundance of sur- 
viving subhalos selected instead by the highest circular 

max,p* 

be- 



vclocity they reached throughout their lifetime, V r 



Subhalos will reach their V r 



max.p 



at a redshift z u 



fore falling into Via Lactea: this type of circular veloc- 



ity selection is designed then to remove the bias intro- 
duced by tidal mass losses, and to highlight the subhalos 
that may have started shining before being accreted by 
their host. Within r2oo, the number of massive galactic 
subhalos that reached a peak circular velocity in excess 
of 10 kms^ 1 at some point during their history is 510, 
about five times larger than their present-day abundance. 
This ratio increases with increasing Vma X ,p and decreas- 
ing radius: a) above a virial temperature T vu - = 10, 000 
K, or a circular velocity V max = 16.7 kms -1 (/i = 0.59 
for fully ionized primordial gas), gas can cool efficiently 
and fragment via excitation of hydrogen Lya. The num- 
ber of subhalos within r2oo that reached this "atomic 
cooling" mass at some point in the past is 135, nearly 
six times larger than their present-day abundance; b) 
within the inner 50 kpc there is today only one subhalo 
with Vmax > 16.7 kms -1 , but there arc 16 surviving rem- 
nants that had this peak circular velocity and were more 
massive at earlier times. If substructure mass regulates 
star formation, then for a given mass threshold many 
more subhalos should have been able to build a sizeable 
stellar mass at some point in the past than indicated by 
their present-day abundance. 

It is important then to investigate the consequences 
of a mass (or circular velocity) cut that picked instead 
the top (say) 65 most massive (or largest Vrnax,p) subha- 
los at all epochs as the hosts of t he kno wn Milky Way 
dwarfs. Following iDiemand et al.l (|2007bf ). such or simi- 
lar samples have been term ed "LBA" (for "larges t be- 
fore accretion" s ubhal os) by IStrigari et alJ (|2007af ) and 
ISimon fc Gehal (|2007h . The idea behind this selection 
is to allow star formation only above a relatively large 
constant critical size, a scenario of permanently ineffi- 
cient galaxy formation in all smallest systems, indepen- 
dently of time-varying changes in the environment like 
those triggered, e.g., by reionization. Today's circular 
velocity distribution of our LBA sample is shown in Fig- 



Madau et al. 



7 




FlG. 5. — Cumulative number of Via Lactea subhalos within r2oo {solid curve) as well a s all Milky Way satell it e galaxies within 42 kpc 
(filled squares), as a function of circular velocity. The data points arc from Mateo (1998), Simon & Gcha (2007), Munoz ct al. (2006), and 
IMartin et al.l (120071) . and assume a maximum circular velocity of V m ax = %/3<r l|Klvpin et al.l[T999t ). The short-dashed curve connecting 
the empty squares shows the expected abundance of luminous satellites after correcting for the sky coverage of the SDSS. Dash-dotted 
curve: circular velocity distribution for the 65 largest Knax.p subhalos before accretion (LBA sample). Long-dashed curve: circular velocity 
distribution for the "fossil of reionization" EF sample. This includes the 61 largest (sub)halos at z = 13.6 [V m ax(z = 13.6) > 4 kms ] 
plus the 4 (sub)halos that reach a Vinax.p > 38 kms -1 after the epoch of reionization and are not in the largest 61 at z = 13.6. 



ure O interestingly, this sample includes 12 of the 14 
subhalos above Vmax = 20 kms -1 identified today, and 
26 of the 35 identified above VJ nax = 15 kms -1 , i.e. the 
most massive today and LBA subpopulations basically 
coincide at large values of Vmax- 5 Therefore a solution 
to the substructure problem in which only the largest 50- 
100 Vmax.p subhalos at all epochs were able to form stars 
efficiently would automatically place the luminous Milky 
Ways dwarfs in the most massive subhalos at the present 
epoch. To match the circular velocity function of the 
LBA sample, however, the observed dwarf spheroidals 
(dSphs) must have circular velocity profiles that peak 
at values well in excess of the stellar velocity disper- 
sion (see Fig. [5] and discussion below). Note that the 
cut in Vmax.p instead of V max of the LBA sample re- 
quires star formation to be inhibited in all subhalos with 

5 Note that the same is not true for the top 10 LBA subhalo s 
(IKravtsov et al.ll200l IDiemand et alJ2007bl ; IStrigari et "aLll2007all , 

as the largest V m ax,p systems suffer the largest mass loss and are 
removed from the top ten list of more massive systems at z = 0. 



Vmax,p < 21.9 kms 1 or virial temperature 

r vir = ^"ffi»".p <i 7 ,000K. (6) 
2k B 

4. SUPPRESSING DWARF GALAXY FORMATION 

The two thresholds for efficient star formation given in 
equations |5]) and ([6]) provide the correct total number of 
luminous Milky Way satellites (assumed to be around 60- 
70), not a match to the observed circular velocity func- 
tion. A careful look at Figure [5] suggests two possible 
solutions to the mismatch problem: 

1 . stars in the Milky Way dSphs are deeply embedded 
within their dark matter halos. The halo circu- 
lar velocity profiles peak well beyond the luminous 
radius at speeds significantly higher that expected 
from the stellar line-of-sig ht velocity dispersion , i.e. 
Vmax ~ 3cr as sug gested bv lStoehr et al.l (|2002t l and 
iPenarrubia et all (|2007l ). This scenario would shift 
the data points in Figure [5] by about a factor s/3 
further to the right, making the mass distribution 
of the luminous Milky Way dwarf spheroidals agree 
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TABLE 1 
The Milky Way satellites 



Name 


Heliocentric distance 


a 


po Core radius r c 


References 




(kpc) 


( km s — 1 ) 


(M Q pc- 3 ) 


(PC) 






Segue 1 


23 


- 


- 


19 


(1) 




Sagittarius 


24 


11.4 


0.07 


550 


(2) 




Ursa Major II 


32 


6.7 


1.14 


81 


(3),(4) 




Willman 1 


38 


1.3 


18 


13 


(5),(6) 




Coma Berenices 


11 


4.6 


2.1 


11 


(1 ,(3 




LMC 


49 


20.2 


- 


- 


(2),(7) 




SMC 


58 


27.5 


- 


- 


(2),(8) 




Bootes 


60 


6.6 


0.34 


115 


9 ,(10) 


Bootes II 


60 


- 


- 


16 


(11) 




Ursa Minor 


66 


9.3 


0.36 


200 


(2) 




Sculptor 


79 


6.6 


0.60 


110 


(2) 




Draco 


82 


9.5 


0.46 


180 


(2) 




Sextans 


86 


6.6 


0.06 


335 


(2) 




Carina 


101 


6.8 


0.17 


210 


(2) 




T Tvqp 1\/Ta m~w T 

KJ ± JCL LY±<j> I yJL ± 


106 


7.6 


0.25 


197 


(3) 




Fornax 


138 


10.5 


0.09 


460 


(2) 




Hercules 


138 


5.1 


0.10 


205 


1),(3) 




Canes Venatici II 


151 


1.6 


0.49 


85 


(1 5,(3) 




Leo IV 


158 


3.3 


0.19 


97 


1),(3) 




Canes Venatici I 


221 


7.6 


0.08 


355 


(3) 




Leo II 


233 


6.7 


0.24 


178 


(2),(12),(15) 


Leo I 


255 


9.2 


0.23 


245 


(13), (16) 


Leo T 


417 


7.5 


0.79 


109 


(3),(17) 


Phoenix 


420 


8.9 


0.14 


310 


(2),(14) 


References. 


— (1) Bclokurov et al 


2007; (2) 


Mateo 1998: (3) 


Simon & Gehal 2007 


; (4) 



eTaLI 12001; (51 IWillman et al.l 120061; (6) IMartin et al.l 120071; (7) 



van der Marel et al.l 



is fe Za ritskv 2006; (9) Bcl okurov et al.112009; ( 10) IMunoz et al.H20061 ; f ll)IWalsh"eTaT 



(12) |Bellazzini et al.||2005f (13) )Mateo et al.ll20C)7l ;Tl4l lYoung et al.H2007l ; (15) IColeman et~ 
(16 nBellazzini et al.ll2Tjoa ; (17) llrwin et al.H2007l . 



2002 



2007 



2007 



with the steep mass function of the most massive 
Via Lactea subhalos today. Its main drawback is 
that it boosts the number of massive inner satel- 
lites to values that are difficult to reconcile with 
ACDM halos. There are five Milky Way dwarfs 
within 50 kpc of the Galactic Center with mea- 
sured stellar velocity dispersions a above 4 kms" 1 : 
the SDSS ultra- faint Coma Berenices, Willman 1, 
and Ursa Major II, together with the LMC and 
Sagittarius. Correcting for the sky coverage of the 
SDSS and using the conversion Vmax ~ 3cr would 
bring the total number of inner satellites above 
Vmax ~ 12 kms -1 to twenty or so. By contrast, 
there are only three subhalos in Via Lactea with 
d < 50 kpc and V mayi > 12 kms -1 today (see the 
right panel of Fig. [3}. Note that, even at the reso- 
lution of Via Lactea, subhalos within (say) 50 kpc 
from the center may still be subject to numerical 
ovcrmcrging. We will compare in detail the ra- 
dial distribution of Milky Way's satellites and Via 
Lactea subhalos in the next section. 

if Knax ~ \/3cr instead, then the discrepancy be- 
tween the number of massive substructure expected 
in ACDM and the known companions of the Milky 
Way can be solved if only approximately one out of 
five subhalos with V max > 20 kms" 1 today houses 
a luminous dwarf satellite. Because of dynamical 
friction and tidal mass losses, the large subhalos 
that are today devoid of stars were typically more 
massive in the past. This scenario therefore re- 
quires star formation to be quenched in objects 
that reached peak circular velocities as large as 



35 — 40 kms 1 in the past. 

To shed further light on the star formation suppres- 
sion mechanism that is acting on galactic substructure 
we have plotted in Figure [5] (left panel) all Via Lactea 
subhalos that reached a V max ,p > 10 kms -1 at redshift 
Zmax- Out of a total of 510 subhalos, 318 reached their 
Vmax.p at redshift z max > 2, 180 at z max > 3, and 80 at 
z > 4. The color coding indicates the different intervals 
of Vmax where those subhalos end up today. Clumps that 
were accreted earlier (z max > 3) are more heavily affected 
by tidal stripping and their V maXiP 's drop by more than 
67%. If stars can form continuously and efficiently only 
above a mass or circular velocity threshold, then clearly 
any quenching mechanism must already be in place at 
early times. 

At a distance of 420 kpc and with a surface bright- 
ness of [iy = 26.9 mag arcsec -2 , Leo T has well- 
measured kinematics from both the stars and the gas, 
a stellar velocity dispersion of 7.5 ± 1.6 kms -1 , an es- 
timated mass just short of 10 7 M Q , an absolute mag- 
nitude My = —7.1, and a total ma s s-to-light ratio 
(M/Ly) of 138 ± 71 (llrwin et all [20071 : iSimon fc Gehal 
I2007t iRvan- Weber et al J 120071) . The detection efficiency 
for satellites this luminous over the SDSS survey area 
is nearly unity, so only a handful of similar objects is 
expected over the whole sky. Even at distances well in 
excess of 500 kpc, objects brigh ter than My = — 8 w ould 
have been found by the SDSS (jKoposov et al.ll2007f ). In 
the right panel of Figure [5] we have plotted in the same 
V maXiP — £ max plane (but on a different Vmax.p scale) the 
sample of 575 halos that lie today in the "field" between 
389 kpc (^200) and 584 kpc (1.5r2oo) an d are more mas- 
sive than 10 7 Mq. And while this numerous subpopula- 
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Fig. 



Le/i panel: All (510) Via Lactea subhalos that reached a Vrr 



> 10 kms 1 at rcdshift 



plotted in the V m ax, P - 2max 

plane. The color coding indicates the different intervals of V ma x where those subhalos end up today: Vmax > V m ax,p/2 (blue dots), 
V m ax,p/3 < Vmax < Vmax,p/2 (cyan dots), Vmax < Vmax,p/3 (green dots). Clumps that were accreted earlier (z m ax > 3) are more heavily 
affected by tidal stripping and their Vmax,p's drop by more than 67%. Right panel: Same (but on a different V m ax,p scale) for the sample 
of 575 halos that lie today in the "field" between 389 kpc (r2oo) an d 584 kpc (1.5r2oo) arlcl are more massive than 10 7 Mq. Most of these 
are accreted late (z m ax ^- 1) and suffer only moderate mass loss. 



tion may be shining just under the detection threshold 
as ultra-faint dwarfs, it seems fair to conclude that small 
dark matter clumps must be relatively inefficient at form- 
ing stars even well beyond the virial radius. The figure 
shows that most of these "field" halos are newcomers, 
i.e. they are accreted late (z max 1) and suffer only 
moderate mass loss. Contrary to what is invoked for 
more massive dwarf sphcroidal s sinking within 5 0-100 
kpc from the Galactic Center (|Maver et a.1.1 l2007h . gas 
stripping by ram pressure in a hot low-density Galactic 
corona is unlikley to be responsible for the truncation 
of star formation in this subpopulation. Indeed, while 
many of these field halos are former subhalos, i.e. passed 
through the host at some earlier time, only 15% of them 
had pericenters within 50 kpc. 

5. A TOY MODEL OF REIONIZATION SUPPRESSION 

A popular model for the suppression of gas accretion 
and star formation in small halos is the "reionization 
scenario". In this hypothesis the observed dSphs cor- 
respond to the surviving substructure population that 
accreted a significant amount of gas before the epoch 
of reionization at redshift > 10 and achieved a poten- 
tial well deep enough to form stars: after reionization, 
UV photoionization feedback results in the quenching of 
star formation in all newly forming subhalos below some 
mass threshold. While this general model has been in- 
vesti gated by a number o f authors in the recent literature 
e.g. iBullock et al.ll2000l; Benson et al.ll200il [So"mervillel 



2002t iKravtsov et all 120041; iRicotti fe Gnedinl 120051; 



Moore et al.ll2006t iDiemand et al.ll2007at IStrigari et al l 
2007aUSimon fc Gehall2007f ). there are considerable dif- 
ferences in the specifics of its many variants. 

The Wilkinson Microwave Anisotropy Probe (WMAP) 
3-year data require the universe to be fully reionized 
by redshift z = 11 ± 2.5 (Spergel et al. 2006). The 
reionization history of a given galaxy-forming region ap- 
pears to depend strongly on the mass of its present-day 



halo, with higher mass systems formi ng earlier on av- 
erage and thus b eing reionized earli er (|Weinmann et all 
120071) . Following [Moore et all (|2006[ ). let us assume that 
the region around the Milky Way was reionized around 
z = 12 — 13, i.e. when reionization was well advanced 
in the universe as a whole. Cosmological hydrodynamic 
simulations of early structure formation in a ACDM 
universe have found that a typical H2 fraction in ex- 
cess of 200 times the primordial value is produced af- 
ter the collapse of halos with virial temperatures above 
1000 K. This is large enough to efficiently cool the gas 
and allow it to collapse and form stars within a Hub- 
ble time un less significant heating occurs during this 
phase (e.g. iMachacek et all 1200 ll lYoshida et all 120031 ; 
iKuhlen fc Madaull2005D . A virial temperature of 1200 K 
corresponds to a circular velocity of V max = 4 kms -1 . 
Prior to reionization, at a redshift of 13.6, there are 61 
halos with V max (z = 13.6) > 4 kms -1 in the Via Lactea 
simulation that will have a surviving bound remnant to- 
day within r2oo- At this early epoch their masses range 
between 1.5 x 10 6 M Q and 1.6 x 10 8 Mq, and only one 
object in the sample is above the critical mass thresh- 
old for atomic cooling. Let us assume that the observed 
dSphs are fossils of this first generation of small-mass 
galaxies, a hypothesis consistent with the fact that old 
populations with ages > 10 Gyr are the domina nt stel- 
lar c o mponent of Milky Way dSphs (see e.g. iGrebell 
[20001: iDolphin et aT1[200l . A cosmic UV background 
at z > 10 has been shown to have little impacts on dwarf 
galaxy-sized objects with V max <; 10 kms -1 that have 
already grown to a susbstantial overdensity when the 
ionizing radiation tu rns on and can self-shield against it 
(jDiikstra et al.ll2004[ ). The most massive of these early- 
forming objects (EF sample) may then retain a fraction 
of their gas even after reionization, and may undergo a 
later episode of star formation triggered e.g. by tidal 
shocking during their accretion into the parent halo. By 
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contrast, after rcionization breakthrough, photoioniza- 
tion heating leads to a dramatic increase in the tempera- 
ture of the intergalactic medium, and gas cannot accrete 
and cool to fuel star formation in newly forming dark 
matter halos that are not sufficiently massive. The ex- 
act value of this threshold has been the subject of many 
studies, and depends on the redshift of collapse, the am- 
plitude and spectrum of the ionizing background, and 
the ability o f the system to self-shield against UV ra- 
diation (e.g. lEfstathioulll992t iThoul fc Weinberg! 119961; 
iNavarro fe Steinmetzl Il997t iDiikstra et alj 12004 ). Here 
we assume that after reionization gas infall and cooling 
occurs only onto dwarf galaxy-sized objects that reach 
the threshold Vmax, P > 38 kms" 1 . This brings our EF 
sample to a total of 65 surviving subhalos (as in the LBA 
sample). Their present-day circular velocity function is 
shown in Figure and it appears to reproduce reason- 
ably well the velocity distribution of the known Milky 
Way satellites. 

6. SUBHALO RADIAL DISTRIBUTION AND CENTRAL 
DENSITIES 

The spatial distribution of luminous satellites provides 
another challenge to ACDM substructure models. Fig- 
ure [7] compares the normalized cumulative radial distri- 
butions of the known Milky Way dwarfs to that of the 
oldest (EF) and highest V ma x,p (LBA) dark matter sub- 
halos of Via Lactea. For plotting and comparison pur- 
poses, we have placed Leo T and Phoenix within r2oo- 
The observed satellites are clearly more biased to lie at 
small radii than both the LBA and the EF samples. The 
median distance of the observed satellites is 90 kpc, and 
the satellite distribution becomes even more compact af- 
ter correcting for the SDSS sky coverage (sec Fig. [7]). 
This must be compared to a median distance of 110 kpc 
for the LBA and EF samples, and to 235 kpc for all 
Via Lactea subhalos. About 40% (50% after SDSS cor- 
rection) of all Milky Way dwarfs are found within 60 
kpc from the Galactic center, compared to only 20% for 
the LBA and EF samples, and to 3% for the subhalo 
population as a whole. The largest subhalos today are 
found at significantly larger radii than both the LBA 
and the EF subpopulation. A s imilar discrepancy has 
bee n previously pointed out by iKravtsov et alj (|2004f ) 
and lWillman et all (|2004h . and may only partially be re- 
solved by incompleteness in the census of distant Milky 
Way dwarfs or by numerical overmerging of subhalos in 
the inner 50 kpc of Via Lactea. Note that the luminous 
satellites appear to be biased towards the center even 
compared to the overall Via Lactea's dark matter distri- 
bution. 

The study of the central densities of dark-matter dom- 
inated dSphs offers another test of ACDM galaxy for- 
mation models and of the systematic small-scale prop- 
erties - only minimally modified by baryonic effects - 
of their halos. To estimate the central mass density, 
po, of pressure-supported low luminosity dwarfs, a "core 
fitting" or King's analysis is usu ally performed ([Kind 
119661 iRichstone fc Tremainelll986l ). Under the assump- 
tion of spherical symmetry, close to dynamical equilib- 
rium, mass following light, and isotropic stellar velocity 
dispersion, one has 

166?7(T 2 



where a is the observed central velocity dispersion in 



kms 



the King profile core radius in pc, and 77 



Po 



M pc- 



CT) 



is a numerical parameter very close to unity. In two- 
component isotropic models in which luminous matter is 
negligible compared to dark matter and in which the core 
radius of the dark matter is larger than that of the visi- 
ble material, equation (J7]) can overes timate the true cen- 
tral d ensity by up to a factor of 2.2 (jPrvor fc Kormendvl 
Il990f ). We have evaluated the central densities of Milky 
Way's spheroidals from equation ([7]) (setting ?] = 1) 
using data in the literature. For the ultra-faint new 
dwarfs, King core radii have been estimated from pub- 
lished Plummer (half- light) radii rp using r c = 0.64 rp 
(|Simon fc Gehal |2007| ). The radii and velocity disper- 
sions assumed for this calculation, as well as our de- 
rived central densities, are given in Table 1. Central 
densities range from < 0.1 M pc~ 3 (e.g. Hercules, Sex- 
tans, Canes Venatici I) to ~ 1 M Q pc~ 3 (Ursa Major 
I, Coma Berenices, Leo T) to ~ 20 M Q pc" 3 (Willman 
1). For comparison, we have computed the mean dark 
matter densities within 300 pc (about 3 times our force 
resolution) from their center, pDM; for all Via Lactea's 
subhalos. These are plotted in the left panel of Fig- 
ure [8] as a function of peak circular velocity for subha- 
los at distances R < 50 kpc, 50 < R < 200 kpc, and 
200 < R < r 2 oo- These central densities increase steeply 
with circular velocities. In the smaller systems a distance 
of 300 pc is comparable to the radius where the peak cir- 
cular velocity is reached, i.e. we probe scales where halo 
profiles are nearly isothermal and indeed the small sys- 
tems follow the pdm(< 300pc) oc V^ax scaling expected 
for isothermal halos quite well. The relation flattens for 
larger subhalos since their scale radii are larger than 300 
pc and our measuments now probe the inner regions, 
which arc shallower than isothermal. 

The right panel of Figure [5] shows a comparison be- 
tween the dark matter densities of subhalos in the LBA 
and EF sample and the central densities of dSphs. Via 
Lactea is clearly starting to resolve the inner dense re- 
gions of Milky Way-like satellites. In particular, the cen- 
tral densities of dSphs with measured core radii in excess 
of 250 pc (i.e. comparable to the radii within which pdm 
is measured) appear to be in good agreement with the 
values found in subhalos of comparable circular veloci- 
ties. A few satellites with core radii smaller than 100 
pc, like Ursa Major II, Coma Berenices, and Willman 1, 
have central densities in excess of 1 M pc~ 3 . The pre- 
dicted dark matter densities we obtain by extrapolating 
inward the steep inner profiles of subhalos arc perfectly 
consistent with the observed values. The left panel of 
the same figure also demonstrates that the typical cen- 
tral densities at a given Vmax are systematically larger for 
inner subhalos. This is a consequence of the larger aver- 
age tidal mass loss suffered by this subpopulation: tidal 
stripping can significantly reduce VJ nax without chang- 
ing the central density, and inner subha los end up hav- 
ing mu ch larger typical concentrations (jDiemand et alJ 
l2007bf) . This concentration-radius relation should be 
taken into account when the central densities of lumi- 
nous dwarfs (constrained by stellar kinematics) are used 
to infer the peak circular velocities of the host subhalos. 

A simple way for estimating the range of host sub- 
halo sizes corresponding to a given Milky Way satellite 
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Fig. 7. — The normalized cumulative radial distribution of Milky Way satellites (lower solid curve). The upper solid curve connecting 
the squares shows the expected distribution of luminous satellites after correcting for the sky coverage of the SDSS. For plotting and 
comparison purposes, we have placed Leo T and Phoenix within r200 = 389 kpc. Long-dashed curve: distribution for the 65 largest V ma x,p 
subhalos before accretion (Via Lactea LBA sample). Dot-short dashed curve: "fossil of reionization" EF sample. Short- dashed-long dashed 
curve: 65 largest Vmax Via Lactea subhalos today. Dotted curve: all 4021 Via Lactea subhalos with Af au b/Mh os t > 10~ 6 . Short-dashed 
curve: dark matter distribution. 



is to use the integrated dSph "central mass" within a 
physical radius comparable to the extent of the lumi- 
nous region, and compare it to those of Via Lactea's 
subhalos at similar galactocentric distances. The cen- 
tral masses Mo. 6 (mass within 0.6 kpc) and Mo.i (mass 
within 0.1 kpc) arc robustly constrained by kinematic 
data. Given Via Lactea's force resolution of 90 pc, 
Mo. 6 is also robustly determined in simulated subhalos. 
Mo.i, however, is likely affected by gravitational soften- 
ing and may be systematically underestimated. Here we 
use the most likely values and 90% confidence regions 
of Mq.i from IStrigari et "all (|2007bf ) and of M . 6 from 
IStrigari et al.l (l2007af). as well as the best-fit M .6 values 
from I Walker et all (|2007r i. We then find all Via Lactea 
subhalos within 40% in distance and within 25% of the 
M .i or M .6 values of each satellite, and then average 



over their V max and V m 



to find the most likely host 



size. The range of plausible peak circular velocities was 
determined using all subhalos within the 90% confidence 

6 For the dSphs marked with an asterisk in Table [5] no subhalos 
were found in the quoted range, and an estimate of the most likely 
host size was made using the subhalo with the smallest quadratic 
sum of the relative difference in distance and central mass. 



region in Mo.i (or Mo. 6)- The results of this exercise 
are given in Table [2] Consistent with the expectations 
from the concentration-radius relation, for a given dSph 
central mass the most likely Vmax values of its host sub- 
halos are larger at large galactocentric distances. If, e.g. 
Ursa Minor with Mo.o = 5.3 M Q were to lie at a distance 
of 300 kpc, its most likely host would have V maXiP = 
31 kms -1 and Vmax = 26 kms -1 . Coma Berenices with 
M .i = 1.9 M Q at 300 kpc would typically be hosted by 



a Vm 



30 kms and V m 



21 kms 1 subhalo. 



7. SUMMARY AND CONCLUSIONS 



In this paper we have reported results from the two 
highest resolution simulations of Galactic CDM substruc- 
ture to date, Via Lactea and le8Ell. The two runs fol- 
low the hierarchical assembly of a Milky Way-sized halo 
down to z = and of an elliptical-sized halo down to 
z = 0.5, respectively, using different cosmological param- 
eters. Our main findings can be summarized as follows: 

1. Via Lactea and le8Ell are characterized by simi- 
larly steeply rising subhalo mass functions. In the 
range 200m p < M su b < 0.01Mh os t, the best-fit 
slope of the differential distribution, dN/dM su h oc 
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Fig. 8. — Left panel: Central dark matter densities (measured within the inner 300 pc) of Via Lactea subhalos as a function of their 
peak circular velocity today. Triangles: all subhalos within 50 kpc from the center. Squares: all subhalos between 50 and 200 kpc. Circles: 
all subhalos between 200 kpc and r2oo- Right panel: Derived central densities of Milky Way dSphs, po = 166<r 2 /r 2 , where a is the stellar 
velocity dispersion in kms~ 1 and r c is the King profile core radius in pc (see text for details). Diagonal crosses: dSphs with r c < 250 pc. 
Filled traingles: dSphs with r c > 250 pc. Empty circles: EF subhalo sample. Empty squares: LBA subhalo sample. 



TABLE 2 

Likely host subhalos of Milky Way satellites 



Name 


Distance 


a 


Central mass a 




Vmax 


References 




kpc 


km s 


10 7 Mq 


kms -1 


kms -1 




Ursa Major II* 


38 


6.7 


0.31+0f 


50[15-67] 


11[11-21] 


(3) 


Willman 1 


13 


1.3 


13+° 15 


37 [ 9-50] 


14 [ 5-17] 


(3) 


Coma B 


45 


4.6 


0.19+°f 


41 [14-63] 


16 [ 8-18] 


(3) 


Ursa Minor 


66 


9.3 


23+ ' 19 


28 [14-63] 


17 [ 8-21] 


(3) 


Sagittarius* 


16 


11.4 


n 7 +2U 
z '-27 


94 


40 


(1) 


Ursa Minor 


66 


9.3 


r n + 1.3 
-°-1.3 


49 [35-63] 


20 [17-22] 


(1) 


Draco 


79 


9.5 


4 9+1-4 


51 [30-67] 


21 [21-23] 


(1) 


Draco* 


79 


9.5 


6.9 


35 


22 


(2) 


Sextans 


82 


6.6 


9+ ' 4 
u - 3 -0.3 


14 [ 8-23] 


7 [ 6-9] 


(1) 


Sextans 


82 


6.6 


2.5 


25 


14 


(2) 


Sculptor 


86 


6.6 


9 7 +0.4 


26 [21-29] 


15 [12-17] 


(1) 


Sculptor 


86 


6.6 


4.3 


32 


23 


(2) 


Carina 


101 


6.8 


04+0.7 


29 [28-30] 


18 [14-23] 


(1) 


Carina 


101 


6.8 


2.0 


25 


12 


(2) 


Fornax 


138 


10.5 


4 3+ 2 - 7 


42 [21-63] 


23 [15-26] 


(1) 


Leo II 


205 


6.6 


2 1+ 1 ' 6 


21 [9-33] 


13 [8-18] 


(1) 


Leo I 


250 


8.8 


4 3+!-6 


29 [21-35] 


23 [15-26] 


(1) 



Refe rences. — (1) IStrigari etaLl [|2007al h- (2) IWalker et all p007l V. (3) IStrigari et"aTI 
(20073). 

a Equal to Mq.i in the top 4 systems and to Mo. 6 in all others. The quoted ranges are 90% 
confidence level regions. For these dSphs no subhalos were found in the quoted range, and 
an estimate of the most likely host size was made us described in the text. 



M sub> is a = -1.86 ± 0.02 for le8Ell and a = 
— 1.90 ± 0.02 for Via Lactea. In the same mass 
range the cumulative mass function has slope —0.92 
for le8Ell and —0.97 for Via Lactea. Compared 
to Via Lactea, lc8Ell produces nearly a factor 
of two more subhalos with large circular veloci- 
ties. When normalized to the maximum circu- 
lar velocity of the host, the substructure circu- 
lar velocity distribution yields 170 subhalos above 



V max = 22.5 kms" 1 = 0.056 Vh ost in le8Ell, com- 
pared to 110 above V max = 10 kms -1 = 0.056 Vhost 
in Via Lactea; 

2. the fraction of the Via Lactea's host halo mass 
brought in by subhalos that have a surviving bound 
remnant today within V2oo and with peak circular 
velocity V max > 2 kms -1 is 45%. The total mass 
accreted in identifiable subunits exceeds 97%, as 
many major mergers at early times left no surviv- 
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ing z = subhalos. Hence at our high resolution 
most of the Via Lactea mass is acquired in dis- 
crete clumps, with no significant contribution from 
smooth accretion. There are (14, 109, 873) subhalos 
today with V max > (20,10,5) kms" 1 , which con- 
tributed (19,30,39)% of the mass of Via Lactea. 
Only 9-10% of the total material they brought 
in remains in self-bound surviviving lumps at the 
present-epoch, the rest having been removed by 
tidal effects to become part of the "smooth" halo 
component. Mass losses are more severe for larger 
subunits. 

3. because of tidal mass loss, the number of Via 
Lactea subhalos that have a surviving bound rem- 
nant today and reached a peak circular velocity of 
> 10 kms -1 throughout their lifetime exceeds half 
a thousand, five times larger than their present-day 
abundance and more than twenty times larger than 
the current (incomplete) tally of observed dwarf 
satellites of the Milky. The 20 largest V max sub- 
halos today have lost between 30% and 99.5% of 
the maximum mass that had before accretion. If 
substructure mass regulates star formation, then 
for a given mass threshold many more subhalos 
should have been able to build a sizeable stellar 
mass at some point in the past than indicated by 
their present-day abundance. 

4. unless the circular velocity profiles of Galactic 
satellites peak at values significantly higher that 
expected from the stellar line-of-sight velocity dis- 
persion, only one out of five subhalos with V max > 
20 kms -1 today must be housing a luminous Milky 
Way dwarf. Any mechanism that suppresses star 
formation in small halos must start acting early, 
at redshift z > 3. Nearly 600 halos with masses 
greater than 10 7 M Q are found today in the "field" 
between r2oo and 1.5r2oo, i-c small dark matter 
clumps appear to be relatively inefficient at form- 
ing stars even well beyond the virial radius. Even 
if most subhalos appear to have no optically lumi- 
nous counterparts, such a large substructure popu- 
lation may be detectable via flux ratio anomalies in 
stron g gravitational lenses (e.g. iMetcalf fc Madaul 
|2001[ ), through its effects on stellar streams (e.g. 



Bcllazzini, M., Gennari, N., & Ferraro, F. R. 2005, MNRAS, 360, 
185 

Bcllazzini, M., Gennari, N., Ferraro, F. R., & Sollima, A. 2004, 

MNRAS, 354, 708 
Bclokurov, V., et al. 2006, ApJ, 647, Llll 
Belokurov, V., et al. 2007, ApJ, 654, 897 

Benson, A. J., Lacey, C. G., Baugh, C. M., Cole, S., & Frcnk, C. S. 

2002, MNRAS, 333, 156 
Bergstrom, L., Edsjo, J., Gondolo, P., & Ullio, P. 1999, PhRvD, 

59, 043506 
Bertschinger, E. 2001, ApJSS, 137, 1 

Bullock, J. S., Kravtsov, A. V., & Weinberg, D. H. 2000, ApJ, 539, 
517 

Bullock, J. S., et al. 2001, MNRAS, 321, 559 

Coleman, M. G., Jordi, K., Rix, H.-W., Grebel, E. K., & Koch, A. 

2007, AJ, 134, 1938 
Colin, P., Avila-Reese, V., & Valenzuela, O. 2000, ApJ, 542, 622 
Dckcl, A., & Silk, J. 1986, ApJ, 303, 39 
Diemand, J., Kuhlen, M., & Madau, P. 2006, ApJ, 649, 1 
Diemand, J., Kuhlen, M., & Madau, P. 2007a, ApJ, 657, 262 



llbata et all l2002t iMaver et all |2002| ). or possibly 
via 7- rays from dark matter annihilation in th eir 
cores (|Bergstrom et al.lll999l ; [Stoehr et al.ll2003t) ; 

5. the observed Milky Way satellites approximately 
follow the overall dark matter distribution of Via 
Lactea. They appear to be more biased to lie at 
small radii than both the LBA and the EF sam- 
ples. The median distance of the observed satellites 
is 90 kpc, and the satellite distribution becomes 
even more compact after correcting for the SDSS 
sky coverage. This must be compared to a median 
distance of 110 kpc for the LBA and EF samples, 
and to 235 kpc for all Via Lactea subhalos. About 
40% (50% after SDSS correction) of all Milky Way 
dwarfs are found within 60 kpc from the Galactic 
center, compared to only 20% for the LBA and EF 
samples, and to 3% for the subhalo population as 
a whole. Such discrepancies may only be partially 
resolved by incompleteness in the census of distant 
Milky Way dwarfs or by numerical overmerging in 
the inner 50 kpc of Via Lactea. 

6. the study of the central densities of dark-matter 
dominated dSphs offers another test of ACDM 
galaxy formation models and of the systematic 
small-scale properties - only minimally modified by 
baryonic effects - of their halos. Via Lactea subha- 
los have central matter densities that increase with 
Vmax and reach pdm = 0.1 — 0.3 M Q pc~ 3 , compa- 
rable to the central densities inferred in Milky Way 
dwarfs with core radii > 250 pc. 
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